Two-Particle Dispersion in Model Velocity Fields 
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We consider two-particle dispersion in a velocity field, where the relative two-point velocity scales 
according to u^(r) oc r" and the corresponding correlation time scales as r(r-) oc r^ , and fix a = 2/3, 
as typical for turbulent fiows. We show that two generic types of dispersion behavior arize: For 
a/2-|-/3 < 1 the correlations in relative velocities decouple and the diffusion approximation holds. In 
the opposite case, a/2 -I- /3 > 1, the relative motion is strongly correlated. The case of Kolmogorov 
fiows corresponds to a marginal, nongeneric situation. 

PACS No: 05.40.-a; 47.27.Qb 



Since the seminal work of Sir L.F.Richardson on par- 
ticles' dispersion in atmospheric turbulence [I] a large 
amount of work has been done in order to understand 
the fundamentals of this process (see [2] and [3] for re- 
views). Based on empirical evidence, Richardson found 
out that the mean square distance i?^(i) = (r^(t)) be- 
tween two particles dispersed by a turbulent flow grows 
proportionally to . The works of Obukhov and Batch- 
elor have shown that Richardson's law is closely related 
to the Kolmogorov-Obukhov scaling of the relative ve- 
locities in turbulent flows. Scaling arguments based on 
dimensional analysis allow then to understand the over- 
all type of the behavior of B?{t), but a full theoretical 
picture of the dispersion process is still lacking [3] . 

The theoretical description of dispersion processes 
typically starts from models, in which one fixes the 
spacial statistics of the well developed turbulent flow 
(Kolmogorov-Obukhov energy spectrum), and discusses 
different types of temporal behavior for the flows [3]. 
Three situations have been considered in-detail so far. 
Here, the white-in time flows represent a toy model which 
allows for deep analytical insights [4-6]. In connection 
with "real" turbulence two other cases are widely dis- 
cussed. One of them supposes that the temporal decor- 
relation of the particles' relative motion happens because 
the pair as a whole is moving with a mean velocity rela- 
tive to an essentially frozen flow structure (as proposed 
by a Taylor hypothesis) [7,8]. Another premise con- 
nects this decorrelation with the death and birth of flow 
structures (" eddies" ) , whose lifetime is governed by Kol- 
mogorov's universality assumption [8,9]. Both these sit- 
uations are extremely awkward for theoretical analysis. 

In the present letter we address the following question: 
What are the generic types of two-particle dispersion be- 
havior in a velocity field whose statistical spatial struc- 
ture is fixed (and similar to that of a turbulent flow), if 
its temporal correlation properties change. This question 
will be discussed in the framework of numerical simula- 
tions and scaling concepts. As we proceed to show, two 
generic types of behavior arise. Thus, the white-in-time 
flow and the Taylor-type situation belong to the classes 
of diffusive and ballistic behavior, respectively. The case 
of Kolmogorov temporal scaling represents a borderline 



situation. 

Let us consider modes of particles' separation in 
a velocity fleld whose two-time correlation function 
of relative velocities behaves as (v(r, ii)v(r, t2)) oc 
(u^(r)) g [(^2 — ii)/T(r)], where T(r) is the distance- 
dependent correlation time. The g-function is deflned 
so that 5(0) = 1 and g{s)ds = 1. The mean square 
relative velocity and the correlation time scale as 



oc Vn 



and 



r(r) oc To 
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One can visualize such a flow as being built up from 
several structures (plane waves, eddies, etc., see Ref.[3]), 
each of which is characterized by its own spatial scale and 
its scale-dependent correlation time. In well-developed 
turbulent flows one has v'^{r) oc £2/3^2/3^ where e is the 
energy dissipation rate, so that a = 2/3. The white-in- 
time flow corresponds to /? = 0. Kolmogorov scaling im- 
plies (3 = 2/3 and Taylor's frozen-flow assumption leads 
to (3 = 1. 

In our simulations we model the two-particle rel- 
ative motion using the quasi-Lagrangian approach of 
Ref. [9] . Parallel to Ref. [9] we conflne ourselves to a two- 
dimensional case, which is also of high experimental in- 
terest [12, 13]. The relative velocity v(r, t) = V x ri(r, t) 
is given by the quasi-Lagrangian stream function rj. This 
function is built up from the contributions of radial oc- 
taves: 
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where ki = 2', and the flow function for one-octave 
contribution in polar coordinates (r, 9) is given by 
r],{hr,t) = F{hr){A,{t)+B,{t)cos{2e + ((>,)). The ra- 
dial part F{x) obeys F{x) = x^(l — x) for < a; < 1 
and F{x) = otherwise, and 0i are quenched random 
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phases. Moreover, Ai{t) and Bi(t) are independent Gaus- 
sian random processes with dispersions = = 
and with correlation times = 2~^^tq. At each time step 
these processes are generated according to Xi [t + At) = 

^Jl ~\At/Ti)^X^{t) + {At/Ti)vlC, where X is A or B, and 
C is a Gaussian random variable with zero mean and unit 
variance. The values of tq and the integration step At 
are to be chosen in such a way that tn ^ At. Typically, 
values of At ~ 10^'* are used. For the noncorrelated flow 
{(3 = 0) the values of Ai{t) and Bi{t) are renewed at each 
integration step At. In the present simulations TV = 16 
was used. The value vo — 1 was employed in the ma- 
jority of simulations reported here, so that only the use 
of a different vq value will be is explicitly stated in the 
following. 

The values of R^{t) obtained from 3000 realizations of 
the flow for several values of /? e [0, 1] are plotted on 
double logarithmic scales in Fig.l, where tq = 0.15 is 
used. One can clearly see that for all /3 a scaling regime 
R^{t) oc f appears. We note moreover that the curves 
for f3 = 0.67, 0.8, 09 and 1.0 are almost indistinguishable 
within statistical errors. The values of 7 as a function of 
f3 are presented in the insert, together with the theoreti- 
cally predicted forms, vide infra. 




Fig. 1. Mean square displacement B?(t) plotted in double 
logarithmic scales. The four lower curves correspond to /3 = 0, 
0.17, 0.33 and 0.5 (from bottom to top). The dashed lines indi- 
cate the slopes 1.5 and 3. The four upper curves for /3 = 0.67, 
0.8, 0.9 and 1 are hardly distinguishable within the statistical 
errors of the simulations. The insert shows the values of 7(/3) 
determined by a least-squares-fit within the scaling region of 
each curve. The error bar shows the typical accuracy of ^1 
7- values. pJThe full lines give the theoretical predictions, Eq.i 
and Eq.(B 



The regimes of dispersion found in simulations can be 
explained within the framework put forward in Ref.[10]. 



The discussion starts by considering l{r) = v{r)T{r), the 
mean free path of motion at the distance r. If this mean 
free path always stays small compared to r, the rela- 
tive motion exhibits a diffusive behavior with a position- 
dependent diffusion coefficient, K{r) oc P{r)/T{r) oc 

1/2 

j,Q!+/3 Xaking as a scaling assumption r oc (r^(t)) = 
-R, one gets that the mean square separation R grows as 
oc f with 
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On the other hand, if l{r) is of the order of r, the mean 
separation follows from the integration of the ballistic 
equation of motion ^i? = v{R) oc i?"/^, see Ref.[ll]. 
Thus, in a flow where a considerable amount of flow lines 
of relative velocity are open, one gets R^ (x V , with 
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The occurrence of either regime is governed by the value 
of the (local) persistence parameter of the flow. 



Ps{r) — l{r)/r — v{r)T{r)/r. 



(6) 



Small values of Ps correspond to erratic, diffusive mo- 
tion, while large values of Ps imply that the motion is 
strongly persistent. The value of the persistence param- 
eter scales with r as Ps{r) oc 7-"/2+/3-i^ Since under par- 
ticle's dispersion the mean interparticle distance grows 
continuously with time, the value of Ps decreases con- 
tinuously for a/2 + f3 < 1, so that the diffusive approx- 
imation is asymptotically exact. For a/2 + P > 1 the 
lifetimes of the structures grow so fast that the diffusive 
approximation does not hold. This situation is one ob- 
served in our simulations for /3 > 2/3. The strong ballis- 
tic component of motion implies that the velocities stay 
correlated over considerable time intervals. The results of 
Fig.l confirm that j{(3) behaves accordingly to Eq.(^ for 
/3 < 2/3 and Eq.(|) for /? > 2/3. We note here that the 
parameters of the simulations presented in Fig.l (vq = 1, 
Tq = 0.15) were chosen in a way that allows to show all 
curves within the same time- and distance intervals. This 
leads to a somehow restricted scaling range and to slight 
overestimate of 7-values in the diffusive domain. 

Strong differences between the diffusive and the ballis- 
tic regimes can be readily inferred when looking at typi- 
cal trajectories of the motion, such as are plotted in Fig. 2 
for the cases /3 — 0.33 and /3 — 0.67. The difference be- 
tween the trajectories is evident both in the (a;,y)-plots 
and in the r(t)-dependences. The curves for f3 — 0.33 
exhibit a random- walk-like, erratic behavior, while the 
curves for /3 = 0.67 show long periods of laminar, di- 
rected motion. In order to quantitatively characterize 
the strength of the velocity correlations we calculate the 
backwards-in-time correlation function (BCF) of the ra- 
dial velocities, as introduced in Rcf.[12]. This function is 
defined as Crir) — {vr{t ~ T)vr{t)) /(w^(t)) and shows, 



2 



what part of its history is remembered by a particle in 
motion. The function is plotted in Fig. 3 against the di- 
mensionless parameter d = —r/t. The functions (ob- 
tained in 10^ realizations each) are plotted for 4 different 
sets of parameters. Here the dashed lines correspond to 
(3 = 0.33, in the diffusive range, for t = 10"^^ 3 • 10"^^ 
10^^ and 3-10^^. These BCF do not scale and are rather 
sharply peaked close to zero, thus indicating the loss of 
memory. The two sets of full lines indicate Cr (r) in Kol- 
mogorov flows, for t = 10~^, 3-10~^, and 10~^. The lower 
set corresponds to the value tq — 0.05 and the upper set 
to the value tq = 0.15. In both cases the functions show 
scaling behavior. No considerable changes in the BCF's 
form occur when further increasing the value of tq up to 
To = 1, thus indicating that the data tq = 0.15 corre- 
spond already to a strongly correlated regime. The form 
of these curves resembles closely the experimental find- 
ings of Ref. [12]. The BCFs for (3 = 1.0 show an overaU 
behavior very similar to the one in Kolmogorov's case. 
Note that as the time grows the curves for (3 — l.Q ap- 
proach those for (3 = 2/3 and probably tend to the same 
limit. The curves for /? = 1.0 and tq = 1 (not shown) fall 
together with those in Kolmogorov's case with tq = 0.15. 
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Fig. 2. Typical trajectories in the diffusive regime (/3 = 0.33) 
and in the Kolmogorov regime (/3 = 0.67). The upper pictures 
show the trajectories in the (x, j/)-plane, the lower ones repre- 
sent the corresponding r(f)-behavior. Note that the scales of 
the right and of the left graphs differ by a factor of 5. 



The similarity in the properties of dispersion processes 
in a Kolmogorov situation with larger Ps (larger tq) and 
in ballistic regime can be explained based on the behav- 
ior of the effective persistence parameter. In the diffu- 
sive regime we supposed that the correlation time of the 
particles' relative velocity scales in the same way as the 
Eulerian lifetime of the corresponding structures. On 
the other hand, in the ballistic regime, (3 > 1 ~ a/2, 
the lifetimes of the structures grow so fast that no con- 
siderable decorrelation takes place during the time the 
particles sweep through the structure. The Lagrangian 



decorrelation process is then connected not to Eulerian 
decorrelation, but to sweeping along open flow lines. 
The effective correlation time then scales according to 
Ts{r) oc r/v{r) cx and the effective value of /3 

stagnates sA (3 — \ — a/2. Thus, all long-time correlated 
cases belong to the same universality class of strongly- 
correlated flows, as the Kolmogorov flows with large Ps, 
for which Eq.(||) and (||) coincide. 
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Fig. 3. The BCF of relative velocities as a function of the 
dimensionless time lag —r/t. The lower group of dashed 
lines corresponds to /3 = 0.33 (the values of t are 0.01, 0.03, 
0.1 and 0.3, from top to bottom). The two groups of full 
curves corresponds to the Kolmogorov case (three curves for 
t = 0.01,0.03,0.1 each). The dotted curves correspond to 
(3 = 1 for the same values of time, see text for details. 



For Kolmogorov flows the ballistic and the diffu- 
sive mechanisms lead to the same functional form of 
i?^(i)-dependence. The functional form of the depen- 
dence of R^{t) on parameters of the flow is oc 



[vQTo/rQ'^^j V in the diffusive situation {Ps <^ 1) and 

i?^ (X (vq/tq^'^^ f in the balhstic case {Ps » 1). As- 
suming that Ps is the single relevant parameter gov- 
erning the dispersion we are lead to the form R^{t) oc 

f{Ps) (^vq/vq^"^^ V , where f{Ps) is a universal function 

of Ps, which behaves as Ps^ for Ps ^ 1 and tends to a 
constant for Ps ^ 1. Thus, for a flxed spatial structure 
of the flow, the following scaling assumption is supposed 
to hold: 



{votr 



= F{vqTq), 



(7) 



which scaling can be checked in our case by plotting 
R^{t)/{vQtY against vqTq. The corresponding plot is 
given in Fig. 4, where we flx i = 0.1, and plot the results 
in three series of simulations. Each point corresponds to 
an average over 5-10^ runs. Here the squares correspond 
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to ?;o = 1 and to the values of tq ranging between 0.01 
and 0.15, the triangles correspond to vq = 0.3 and th tq 
between 0.033 and 0.5, and the circles to tq = 0.1 and to 
values of Wo between 0.1 and 1.5. The error bar indicates 
a typical statistical error as inferred from 5 similar series 
of 5 • 10* runs each. The scaling proposed by Eq.(0) is 
well-obeyed by the results. Some points outside of the 
range of Fig. 4 were also checked. Thus, for larger val- 
ues of vqTo the values of R"^ {t) / {vot)^ seem to stagnate. 
On the other hand, increasing vqTq to values larger than 
0.3 (i.e. approaching the frozen flow regime) leads to a 
strong increase in fluctuations, making the results less 
reliable. 
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Fig. 4. The values oi R'^{t)/{vot)'^ plotted against votq. The 
dashed line is drawn as a guide to the eye. 



Let us summarize our findings. Thus, we considered 
two-particle dispersion in a velocity field scaling accord- 
ing to u^(r) oc r^^^ and r(r) oc . We show that two 
generic types of behavior are possible: For a/2 + P < 1 
the diffusion approximation holds and the increase in 
the interparticle distances is governed by the distance- 
dependent diffusion coefficient K{r) oc r"^/^ . In the 
opposite case a/2 + P > 1 the relative velocities stay 
strongly correlated. The transition between the two 
regimes takes place exactly for the Kolmogorov flow, for 
which a/2 + (3 — 1. In this case the properties of the 
dispersion process depend on the persistence parameter 
of the flow. 
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